On the Scaling Behavior of the Abelian Sandpile Model 
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The abelian sandpile model in two dimensions does not show the type of critical behavior familar 
from equilibrium systems. Rather, the properties of the stationary state follow from the condition 
that an avalanche started at a distance r from the system boundary has a probability proportional 
to \j\fr to reach the boundary. As a consequence, the scaling behavior of the model can be 
obtained from evaluating dissipative avalanches alone, allowing not only to determine the values of 
all exponents, but showing also the breakdown of finite-size scaling. 
PACS numbers: 64.60.Lx 
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Since its introduction in 1987, the sandpile model has 
been considered as the prototype of a self-organized criti- 
cal (SOC) system jjj. Computer simulations suggest that 
irrespective of the initial conditions and of details of the 
model rules, the system self organizes into a "critical" 
state with a power-law size distribution of avalanches. 
The concept of SOC is thought to explain the frequent 
occurrence of power laws in nature. 

Considerable effort has been made to determine the 
scaling behavior of the two-dimensional abelian sandpile 
model. The advantage of this model is that several of 
its properties can be calculated analytically j2U|]. How- 
ever, the exponent r r characterizing the distribution of 
avalanche radii, and related exponents for the distribu- 
tion of the avalanche duration t and the number of top- 
plings s during an avalanche have resisted any attempt 
of an analytical calculation. The numerical determina- 
tion of these exponents is hampered by the fact that the 
double logarithmic plots show slopes that increase with 
increasing system size, until finite-size effects set in, indi- 
cating that the asymptotic scaling behavior does not yet 
occur for the system sizes accessible to computer simu- 
lations. Thus, the predictions for the value of t s vary 
between 1.22 {§ and 1.27 § or even 1.29 @. The latter 
two results where obtained under the assumption that 
the system displays finite-size scaling (FSS). 

Recently, evidence was found that the abelian sand- 
pile model does not display FSS M. This result was 
obtained from an investigation of multifractal spectra. 
Indeed, there is no a priori reason why the abelian sand- 
pile model should display FSS. The concept of FSS has 
its origin in equilibrium critical phenomena where a small 
finite system cannot be distinguished from a small part 
of a large system. However, this is not the case for the 
abelian sandpile model in two dimensions. In a small sys- 
tem, the sites that participate in an avalanche may topple 
a few times during the duration of the avalanche. In a 
larger system, the number of topplings per site during 
a large avalanche is larger. Thus, locally collected data 
(i.e. the number of topplings of a given site) contains in- 
formation about the size of the system. Also, finite-size 
scaling is based on the assumption that boundaries play 
no special role in the system. However, the boundaries of 
the abelian sandpile model play an essential role as they 



are the only place where sand can leave the system. 

It is the purpose of this letter to elucidate the scaling 
behavior of the abelian sandpile model. Since bound- 
aries play a special role in the system, it is essential to 
consider avalanches that reach the boundaries of the sys- 
tem separately from those that don't. It turns out that 
the dissipative avalanches display beautiful power laws. 
These power laws imply that each avalanche started at 
a distance r from the system boundary has a probability 
proportional to to reach the boundary and, if it 

does so, dissipates on an average y/r sand grains. From 
the scaling behavior of dissipative avalanches, we obtain 
the values of the critical exponents characterizing the sys- 
tem, and we see compelling evidence for the violation of 
FSS. 

The two-dimensional abelian sandpile model is defined 
on a square lattice with L 2 sites. At each site i an inte- 
ger variable z% represents the number of grains. Grains 
are added individually to randomly chosen sites of the 
system. When the number of grains at a site i exceeds 
z c = 3, site i is unstable and topples, its height being 
reduced to Zi — 4, and the heights {zj} of all four nearest 
neighbors being increased by one. If i is a boundary site 
with I < 4 neighbors, 4 — I grains leave the system. If 
a neighbor j becomes unstable due to the addition of a 
grain, it also topples, and the avalanche stops when a new 
stable configuration is reached. During an avalanche, no 
new grains are added to the system. The size s of an 
avalanche is defined to be the total number of topplings. 
The radius of an avalanche is in this paper taken to be the 
maximum distance of toppled sites from the starting site 
of the avalanche. It has proven useful to decompose an 
avalanche into "waves of toppling" ||. The nth wave of 
toppling begins when the starting point of an avalanche 
topples for the nth time, and all those sites belong to 
it that topple immediately after a nearest neighbor that 
belongs to the same wave has toppled. 

Let us first discuss the radius distribution of 
avalanches, n r (r, L). Since its asymptotic (i.e., L — > oo) 
behavior is difficult to extract from simulation data, re- 
cent results for the exponent r r vary between 7/5 R] 
and 5/3 ||. In ||], a version of the abelian sandpile 
model was studied where sand grains were only dropped 
at the center of the system. It was found that the frac- 
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tion \j\fL of all avalanches reach the boundaries. This 
finding would imply that r r = 3/2, just as predicted in 
This is indeed the correct exponent, as a separate 
evaluation of dissipative and nondissipative avalanches 
shows. We denote by rir (r, L) the radius distribution 
of nondissipative avalanches, i.e., those avalanches that 
do not reach the boundaries. The distribution of dissi- 

(2) 

pative avalanches is denoted n r (r, L). Normalization of 
the radius distribution requires 

/oo 
[n r ^{r,L) + n^(r,L)}dr = 1. 

Fig. |l| shows the result of computer simulations for 
n?\r,L). 
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FIG. 1. The radius distribution of avalanches that touch 
the boundaries of the system, for L =128, 256, 512, 1024, 
2048. The thick solid line is a power law with an exponent 
-1/2. 

(2) 

The distribution rir (r, L) has the form 

n^(r : L)^L-^g(r/L). (1) 

(2) 

The prefactor is such that the total weig ht of n r '{r,L) 
is J L rir (r, L)dr ~ l/\fL, in agreement with the result 
in that the fraction 1 / \[~L of all avalanches reach the 
boundaries. Since grains are dropped randomly into the 
system, the probability that an avalanche is triggered at 
a fixed distance r <C L from the boundary is proportional 
to 1/L, leading to n r 2) (r, L) = (l/L)/(r) for r < L. To- 
gether with Eq. ([!]), this gives g(r/L) oc (r/L)^ 1 / 2 for 
small r/L. From Fig. [ij it can be seen that this l/xA 7- 
behavior extends almost over the entire range of r values. 
An avalanche triggered at a distance r from the bound- 
ary consequently has a probability proportional to 1/y/r 
to reach the boundary, and if it reaches the boundary it 
dissipates on an average of the order y/r sand grains. 

Fig. H shows the radius distribution of avalanches that 
do not reach the boundaries. For the system sizes stud- 
ied in the simulations, no power law is visible. The slope 
seems to become steeper with increasing r. (The last part 
of each curve, which begins at the point where it splits 
from the curves for larger system sizes, is due to finite-size 
effects and it must be ignored in the subsequent discus- 
sion which refers to the thermodynamic limit L — > oo.) 

However, the curves rir {f, L) must approach a power 
law rir(r, oo) ~ r~ Tr for sufficiently large r. The reason 
is that the last wave of topplings of an avalanche has been 



proved to be distributed according to a power law r 7 ^ 4 
ifjof , and rir (r, oo) cannot be steeper than this. From 

(2) 

the above results for n r (r, L) , we can even deduce the 
value of 7>: The probability that an avalanche triggered 
at distance r from the boundary reaches the boundary 

(21 

is proportional to Ln r (r,L), and it is identical to the 
probability that an avalanche triggered in the interior of 
the system reaches at least a radius r, which is given 
by rir (r', oo)dr'. The reason is that the landscape 
of height values in the abelian sandpile model does not 
show long-range correlations. Rather, correlations decay 
as fast as r~ 4 ^ , implying that avalanches spread in the 
same way everywhere in the system, as long as they do 
not encounter the boundaries. Consequently, n^\r, oo) 
must fall off as r -3 / 2 for sufficiently large r, fixing the 
value of T r = 3/2. Fig. || shows that this asymptotic value 
is barely reached for the system sizes accessible to com- 
puter simulations. This explains why evaluations based 
on the complete radius distribution ni (r, L) + n r 2 ^ (r, L) 
do not reveal the correct exponent. 
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FIG. 2. The radius distribution of avalanches that do not 
touch the boundaries of the system, for L =128, 256, 512, 
1024, 2048 from left to right. The dotted line is a power law 
with the exponent —3/2. 



The distribution of the number of topplings n s {s, L) = 
TVs (s, L) + n s 2 \s, L) in avalanches can be analyzed in a 
similar way. The superscript (1) refers again to nondis- 
sipative, the superscript (2) to dissipative avalanches. If 
the number of topplings s in an avalanche of radius r was 
always of the order r D , simple scaling would hold, result- 
ing in a FSS form n s (s, L) ~ s~ Ts C(s/ L D ) and a scaling 
relation D = (r r — l)/(r s — 1). The breakdown of simple 
scaling was pointed out for the first time in HQ . It can 
best be visualized by looking at the distribution of the 
number of topplings ni spanning ' ) (s, L) of those avalanches 
that span the entire system, i.e. that touch all four edges 
of the boundary. These avalanches have a radius of the 
order L. 
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FIG. 3. The distribution ni spanning) (s, L) of the num- 
ber of topplings in system-spanning avalanches for L = 
2048,1024,512,256,128 (from widest to narrowest curve). 

Evaluation of the curves represented in Fig. || shows 
that the maximum scales as L D with D = 9/4, the lower 
cutoff as L 2 , and the upper cutoff as L 11 / 4 . These expo- 
nents are obtained by scaling s in such a way that the 
maxima, the left tails, or the right tails collapse. The 
mean scales roughly as L b l 2 . Clearly, a complete scal- 
ing collapse would only be possible if the four quantities 
were characterized by the same exponent. The scaling 
behavior of the mean implies 

/oc 
sn s (s, L)ds ~ L 2 , 

in agreement with the analytical result given in . 

Fig. ^ shows the distributions (s, L) of the number 
of topplings in dissipative avalanches. 




FIG. 4. The distribution n a (s,L) of the number of top- 
plings in dissipative avalanches for L = 
2048,1024,512,256,128,64 (from widest to narrowest curve). 
The dotted line is a power law with the exponent -7/9. 

The curves have the form nk\s) ~ L~ 1 s~ 7 / 9 for 
sufficiently small s, and then have a cutoff that does 
not display FSS but is similar to the one for spanning 
avalanches. The total weight of dissipative avalanches, 

£9/4 y<2\ 

1/yL, is proportional to n s '(s,L)ds. The expo- 



nent t s can be derived from the condition 

/>00 

a x - x t D Ln!?>{a t L)~ J n^(s',oo)ds' , 

giving t s = 11/9, which agrees very well with the value 
obtained by Manna [Q. The asymptotic slope — t s is 
barely visible in nh\s,L) for the system sizes accessible 
to computer simulations. For the value of the exponent 
t s , the multiscaling features are irrelevant, since only the 
exponent D = 9/4 enters the evaluation of t s . However, 
the multiscaling behavior is relevant for the result s = L 2 , 
which stands in no relation to the value of r a . 

The failure of FSS in the system is due to a broad 
distribution of the number of waves of toppling in 
avalanches. Single waves of toppling display FSS, as 
shown in Fig. |^ for the first wave of toppling. 




FIG. 5. The distribution n s (s,L) of the number of top- 
plings in dissipative first waves for L = 2048,1024,512,256,128. 
The dotted line is a power law with the exponent -7/8. 

The scaling with L 2 follows directly from the fact that 
waves are compact and that each site topples once in a 
wave. Repeating the evaluation performed already twice, 
we find an exponent t s = 11/8 for the first wave, which 
is identical to the analytically derived exponent for the 
last wave [jUTf , and different from the exponent 5/4 for 
boundary avalanches in a sector of 360° ||. From the 
condition that waves are compact, it follows immediately 
that T r = 7/4 for the first wave, which means that a frac- 
tion proportional to 1/L 3 / 4 of all first waves reach the 
boundary. The result r r = 7/4 follows also directly from 
a scaling collapse of the radius distribution of dissipative 
first waves. 

For the distribution of duration times t of avalanches 
exactly the same analysis can be performed as for the 
number of topplings s. One finds again a violation of 
FSS, the typical duration time of an avalanche of radius 
r being ~ r z with z = 4/3, the upper cutoff of the dura- 
tion time being larger than r 3 / 2 , and the mean duration 
time of avalanches, averaged over all avalanche sizes, be- 
ing t ~ L. 

The following table summarizes the results for the ex- 
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poncnts: 











z 


D 


I s * wave 
avalanche 


7/4 
3/2 


8/5 
11/8 


11/8 
11/9 


5/4 
4/3 


2 

9/4 



In summary, we have shown that the scaling behavior 
of the abelian sandpile model in two dimensions is tied 
to the condition that the fraction \ j\fr of all avalanches 
started at a distance r from the boundary reach the 
boundary, rather than to a self-similarity of large and 
small avalanches. As a consequence, the scaling proper- 
ties are seen most clearly from the dissipative avalanches, 
and a violation of FSS can occur. 

The condition that a fraction proportional to 1/VL of 
all avalanches reach the boundaries, together with the 
two conditions s ~ L 2 and t ~ L seem so simple that 
they should be derivable from simple argument. In fact, 
s ~ L 2 follows from the diffusive motion of sandgrains 
j| , and appears to be also valid in higher dimensions jll) . 
The first condition is equivalent to the statement that the 
probability of triggering a system spanning avalanche, p, 
is proportional to the density of surface sites, p, that top- 
ple during a system-spanning avalanche. The reason is 
that the product ppL d ~ x is the mean number of grains 
dissipated during an avalanche, which must be identi- 
cal to 1, resulting in p ~ p ~ 

y L (d-i)/2 ( d denotes the 
dimension of the system) . A condition p ~ p holds for ex- 
ample for percolation clusters in systems above the perco- 
lation threshold. It also holds at the percolation thresh- 
old, if there are only a finite number of spanning clusters 
in the system. In three dimensions, this condition would 
lead to r r = 2, which agrees with the simulation result 
given in |^. If avalanches are compact, as suggested by 
computer simulations, the exponent r Q characterizing the 
area distribution in three dimensions is then t q = 4/3, 
giving a mean avalanche area ~ L 2 . Together with the 
condition s — L 2 1 it follows that the mean number of top- 
plings is of the same order as the mean area, leading to 
t s =4/3 and to the conclusion that the number of waves 
of topplings remains finite in the thermodynamic limit, 
and that FSS is not violated in three dimensions. Finally, 



a compact avalanche of s = L 2 sites has a radius L 2 l d . In 
two dimensions, this corresponds to the mean avalanche 
duration t. If we assume the same in three dimensions, 
we have t ~ L 2 ^ 3 , and we find r f = 8/5, again in agree- 
ment with simulation results |6|. In dimensions above the 
upper critical dimension 4 |12|, avalanches are no more 
compact, and the system may contain many spanning 
clusters, and the above arguments can therefore not be 
applied. 

It thus seems that the scaling behavior of the abelian 
sandpile model follows from a few simple principles, a 
connection which still has to be explored in greater depth. 
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